arXiv: 1508.03115vl [astro-ph.HE] 13 Aug 2015 


Astron. Nachr. / AN 999, No. 88, 789 -17931(2006) / DOI please set DOl! 


Ultraluminous X-ray pulsar: accreting magnetar? 
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Magnetars are a special kind of neutron stars. There may also be accreting magnetars. From the studies of isolated mag- 
netars, it is known that a neutron star with a strong dipole field only is not a magnetar. Super-slow X-ray pulsars may 
just be accreting high magnetic field neutron stars. The ultra-luminous X-ray pulsar NuSTAR J095551+6940.8 may be an 
accreting magnetar. It may be an accreting low magnetic field magnetar with multipole field of 10 14 G and dipole field of 
10 12 G. This point of view is consistent with the study of isolated magnetars. An ultra-luminous X-ray pulsar phase in the 
binary evolution may result in massive millisecond pulsars. 
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1 Introduction 

Pulsars are rotating neutron stars. They were first discovered 
in 1967 (Hewish et al. 1986). With the advance of X-ray 
astronomy, accreting neutron stars were discovered in the 
1970s (Giacconi et al. 1971). They manifest themselves as 
accreting X-ray pulsars. In 1982, millisecond pulsars were 
discovered (Backer et al. 1982). They are commonly as¬ 
sumed to be recycled neutron stars via low mass X-ray bi¬ 
naries (Alpar et al. 1982). In this way, normal pulsars, ac¬ 
creting neutron stars, and millisecond pulsars are linked to¬ 
gether. In the 1990s, the idea of strongly magnetised neutron 
stars (i.e. magnetar) were theoretically proposed and obser- 
vationally discovered (Duncan & Thompson 1992; Kouve- 
liotou et al. 1998). Since magnetars are just a special kind 
of isolated neutron stars, it is natural to ask “where are ac¬ 
creting magnetars?” (Woods & Thompson 2006). 

According to the general results of accretion power, the 
accreted matter onto magnetars would release their gravita¬ 
tional energy mainly in X-rays. The typical photon energy 
would be in the range of 1 keV to 50 MeV (Frank et al. 
2002). Therefore, accreting magnetars are expected to be 
some kind of X-ray pulsars. Currently, four kind of X-ray 
pulsars are known. 

1. Rotation-power X-ray pulsars. The X-ray emissions of 
the Crab pulsar and the Vela pulsar etc are assumed to be 
originated from their rotational energy (Becker 2009). 

2. Accretion powered X-ray pulsars. They are accreting 
neutron stars in X-ray binaries (Bhattacharya & van den 
Heuvel 1991). 

3. Magnetars. Observationally, anomalous X-ray pulsars 
and soft gamma-ray repeaters are thought to be mag¬ 
netar candidates (Olausen & Kaspi 2014; Mereghetti et 

* Corresponding author: e-mail: tonghao@xao.ac.cn 


al. 2015; Tong & Xu 2014). The persistent and burst 
energy of anomalous X-ray pulsars and soft gamma-ray 
repeaters may originate from the neutron star’s magnetic 
energy. 

4. If the above three kinds of energy reservoir are all out of 
reach, then the central neutron star can still emit X-rays 
due to its relic thermal energy (Page et al. 2006). The 
thermal X-ray emission may also present in other types 
of X-ray pulsars, e.g. rotation-powered X-ray pulsar. X- 
ray dim isolated neutron stars (Turolla 2009) and central 
compact objects (Gotthelf et al. 2013) have very nice 
thermal X-ray spectra. 

Then both theory and observations should tell the differ¬ 
ences between accreting magnetars and the four kinds of 
X-ray pulsars known at present. 

For accreting compact objects, including accreting mag¬ 
netars, the maximum luminosity for steady spherical accre¬ 
tion is (i.e. Eddington luminosity, Frank et al. 2002): 

^Edd = 1-3 x 10 38 Mi ergs -1 , (1) 

where Mi is the mass of the central object in units of so¬ 
lar masses. For a neutron star with whose radius is about 
10 6 cm, its effective temperature is about 2 x 10' K when 
radiating at the Eddington luminosity. If the central neutron 
star is radiating at hundreds or thousands of times super- 
Eddington luminosity, then its radiation may mainly in the 
hard X-ray range. This is why the giant flares of magne¬ 
tars and ultra-luminous X-ray pulsar are mainly observed 
using hard X-ray telescopes (Mereghetti 2008 and refer¬ 
ences therein; Bachetti et al. 2014). For the ultra-luminous 
X-ray pulsar, it can also be observed in the soft X-ray range. 
However, its timing properties can only be studied using the 
NuSTAR hard X-ray telescope (Bachetti et al. 2014). 

Some basic knowledge of isolated magnetars is pro¬ 
vided in Section 2. The case of super-slow X-ray pulsars 
is discussed in Section 3. The discovery of ultraluminous 
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X-ray pulsar and its explanations in the accreting magnetar 
model are presented in Section 4. Discussion on current ob¬ 
servations and models is given in Section 5. Conclusion of 
this paper is provided in Section 6. 

2 Basic physics of isolated magnetars 

The period and period-derivative diagram of various pulsar¬ 
like objects is shown in figure |T] From the measured period 
and period derivative, the characteristic magnetic field and 
characteristic age of pulsars and magnetars can be obtained 
by assuming magnetic dipole braking in vacuum (Tong 
2015a). For magnetars with typical period of 10 seconds, 
and period derivative of 10 -11 , their characteristic magnetic 
field is 


B c = 3.2 x 10 19 y/pP ~ 3 x 10 14 G, 
and characteristic age 

(2) 

P 

t c = —- ~ 2 x 10 yr. 

2 P 

(3) 


This forms the basic belief that magnetars are young neu¬ 
tron stars with strong magnetic field (Mereghetti 2008). 
However, this basic picture is too simple to be true. Some 
pulsars can have magnetic field as high as 10 14 G while not 
showing significant magnetar activities (Ng & Kaspi 2011). 
The discovery of low magnetic field magnetars (with char¬ 
acteristic magnetic field less than 7.5 x 10 12 G, Rea et al. 
2010) further shows that the dipole magnetic field is not 
the key ingredient of magnetars. Magnetars may be neutron 
stars with strong multipole field only (Tong et al. 2013 and 
references therein). The strong multipole field is responsible 
for the bursts (including giant flares), super-Eddington lu¬ 
minosity during bursts, persistent emissions, and spin down 
of magnetars etc. 

- Giant flares Three giant flares are observed from three 
different magnetars (Mereghetti 2008 and references 
therein). The giant flare is made up of a spike and a 
pulsating tail (lasts for hundreds of seconds). During 
the pulsating tail, the X-ray luminosity can be as high 
as 10 42 ergs _1 which is about 10 4 times the Edding¬ 
ton luminosity. In the magnetar model, the strong mul¬ 
tipole field as high as 10 15 G provides the energy bud¬ 
get of giant flare (Thompson & Duncan 1995). For a 
B ~ 10 15 G multipole field in the vicinity of the neu¬ 
tron star surface, its typical spatial scale may be as that 
of the neutron star radius about R ~ 10 6 cm. Then the 
total magnetic energy will be 

B 2 

£ b ~t-R 3 ~ 4 x 10 46 erg. (4) 

o7T 

If a significant amount of this energy can be release 
is a catastrophic way, it may corresponds to the giant 
flare of magnetars (Yu & Huang 2013). At the same 
time, if the magnetic energy can be release gradually 
this may correspond to its persistent state X-ray lumi¬ 
nosity (Thompson & Duncan 1996). For a typical time 


scale of r ~ 10 5 yr which is the typical age of magne¬ 
tars, the X-ray luminosity is 

~ iO^ergs- 1 . (5) 

T 

Furthermore, the magnetic field will also suppress the 
scattering cross section between electrons and photons 
(Paczynski 1992). The traditional Eddington luminosity 
uses the Thomson cross section and does not consider 
the effect of strong magnetic field (Frank et al. 2002). In 
the presence of strong magnetic field, the corresponding 
critical luminosity is (Paczynski 1992): 


_+_ sa2x (_e_) 4,J 

Lem (id ,2a ) 


(6) 


For a magnetic field strength of 10 15 G, the corre¬ 
sponding critical luminosity can be 10 4 times the tra¬ 
ditional Eddington luminosity. This may explain the 
super-Eddington pulsating tail in magnetar giant flares 
(Paczynski 1992; Mereghetti 2008). If a magnetar is 
accreting from a binary companion, its luminosity can 
also be hundreds or thousands times the traditional Ed¬ 
dington luminosity. This inspired the accreting magne¬ 
tar model for the ultra-luminous X-ray pulsar (Section 4 
below). 

- Timing behaviors Compared with normal rotation- 
powered pulsars, magnetars are always variable (Rea 
& Esposito 2011). They can have burst, and outbursts. 
During these active episode, the magnetar’s timing be¬ 
havior also changes (e.g. Archibald et al. 2015 and ref¬ 
erences therein; Tong & Xu 2014). The characteristic 
magnetic field of the low magnetic field magnetar SGR 
0418+5729 may be related to its small magnetic incli¬ 
nation angle (Tong & Xu 2012). The decreasing period 
derivative of the second low magnetic field magnetar 
Swift J1822.3—1606 may be due to its decreasing par¬ 
ticle wind (Tong & Xu 2013). The anti-correlation be¬ 
tween timing and X-ray luminosity of the galactic cen¬ 
tre magnetar SGR J1745—2900 may due to its chang¬ 
ing geometry (Tong 2015b). The so called anti-glitch of 
magnetar IE 2259+586 may simply be a period of en¬ 
hanced spin down rate (Tong 2014). All these timing 
behaviors of magnetars is understandable in the wind 
braking model of pulsars (Tong et al. 2013). 

- Spectra properties Magnetars can be detected in the 
hard X-ray range (Gotz et al. 2006). There may be a 
cutoff around 1 MeV in the high energy spectrum of 
magnetars (Abdo et al. 2010). A possible cutoff at about 
130 keV is found in magnetar 4U 0142+61 using nine 
years of INTEGRAL data (Wang et al. 2014). A cutoff 
at 130 keV will rule out ultra-relativistic emission mod¬ 
els for the hard X-rays. The putative electrons should be 
at most mildly relativistic. Compared with the hard X- 
ray emission of rotation-powered pulsars and accretion- 
powered pulsars (Mereghetti 2015, figure 3 there; Wang 
et al. 2014), magnetar’s hard X-ray emission will domi¬ 
nate above 100 keV. Therefore, for an accreting magne- 
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tar, a hard X-ray tail above 100 keV may present. This 
may be another signature of accreting magnetars. 

In summary, the study of isolated magnetars tells us 
that the difference between magnetars and normal rotation- 
powered pulsars is their strong multipole field. It is not their 
positions on the period and period diagram which is deter¬ 
mined mainly the magnetospheric torque. Therefore, if peo¬ 
ple want to find accreting magnetars, the signature of strong 
multipole field in accreting systems are expected (Tong & 
Wang 2014): (1) magnetar-like bursts. The magnetar burst 
is different from type I or type II burst in accreting systems 
(Gavriil et al. 2002). (2) a hard X-ray tail above 100 keV. 
China’s future hard X-ray modulation telescope (HXMT) 
may contribute on these aspects (Guo et al. 2015a). 


Their large period derivative also require a very high mag¬ 
netic field. Therefore, super-slow X-ray pulsars are often as¬ 
sumed to be accreting magnetars. However, the magnetic 
field determined from the timing behavior is the large scale 
dipole field and it also depends on the accretion torquePj. 
More importantly, even if the central neutron star has a very 
high dipole field, it is at most an accreting high magnetic 
field neutron star. In order to say that an accreting magne¬ 
tar is seen evidences of strong multipole field are required. 
Considering that super-slow X-ray pulsars may have very 
high dipole field, it is also possible that some of them may 
also have very strong multipole field. Therefore, they are 
candidates for searching of magnetic activities in future ob¬ 
servations. A comparison between isolated magnetars and 
super-slow X-ray pulsars are summarised in Table 1 in Tong 
& Wang (2014). 



Fig.l Period and period-derivative diagram of various 
pulsar-like objects. Black dots are normal rotation-powered 
pulsars, blue squares are magnetars (empty squares are 
radio-loud magnetars). Dot with a circle means the pulsar is 
in a binary system. Updated from figure 1 in Tong & Wang 
(2014). 


3 Super-slow X-ray pulsars are not accreting 
magnetars 

Super-slow X-ray pulsars are a special kind of accreting 
neutron stars (Tong & Wang 2014): (1) They have period 
longer than 10 3 s; (2) large period derivatives 10 -7 — ICC 6 
which may be either positive or negative; (3) their X-ray lu¬ 
minosity (1 — 100 keV) is about 10 34 — 10 37 erg s _1 which 
may be accretion of stellar wind from the neutron star’s high 
mass companion. Theoretically, if their long pulsation pe¬ 
riod is equal to the accretion equilibrium period, the cor¬ 
responding magnetic field required will be about 10 15 G. 


4 Ultraluminous X-ray pulsar and the 
accreting magnetar model 

Bachetti et al. (2014) discovered that NuSTAR 
J095551+6940.8 is an ultraluminous X-ray pulsar with 
isotropic X-ray luminosity about 10 4 °ergs _1 . One previ¬ 
ously identified ulltraluminous X-ray source M82 X-2 may 
be the soft X-ray counterpart of this ultra luminous X-ray 
pulsar (Feng & Kaaret 2007; Kong et al. 2007; Bachetti et 
al. 2014). The pulsation period of the pulsar is 1.37 s. At 
the same time, the pulsar is spinning up at an rate about 
P « —2 x Hr 10 . The super-Eddington luminosity and 
the spin-up trend pose challenges to the accretion theory of 
neutron stars. 

From previous study of isolated magnetars, an aged 
magnetar is more like to be a low magnetic field magnetar 
(Rea et al. 2010; Turolla et al. 2011). For an accreting mag¬ 
netar, the neutron star may be relatively old with age about 
10 6 yr. Then it is an accreting low magnetic field magnetar. 
Tong (2015c) pointed that NuSTAR J095551+6940.8 may 
be an accreting low magnetic field magnetar. The magne¬ 
tar strength multipole field (about 10 14 G) is responsible for 
the super-Eddington X-ray luminosity (cf. equation®). The 
much lower large scale dipole field (e.g. about 10 12 G) is 
responsible for the interaction between the accreting flow 
and the central neutron star. In this case, the spin-up rate 
is understandable even in the traditional form of accretion 
torque. 

Far away from the neutron star, the accretion flow may 
of a disk form if the accreted material posses enough angu¬ 
lar momentum. Closer to the neutron star the effect of the 
central star magnetic field is stronger. The Alfven radius is 
defined as where ram pressure of the accretion flow equals 
the local magnetic energy density (Shapiro & Teukolsky 
1983). Inside the Alfven radius, the accretion flow is con¬ 
trolled by the magnetic field. The Alfven radius depends on 

1 Using a different torque, a strong dipole field may not be required 
(Shakura et al. 2012). 
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the geometry of the accretion flow. It is approximately (Sha- 
pior & Teukolsky 1983; Lai 2014) 

Ra = 3.2 x 10 8 ^ 7 M“ 1 / 7 M“ 2 / 17 cm, (7) 

where /Z 30 is dipole magnetic moment of the central neu¬ 
tron star, Mn is the mass accretion rate in units of 
10 17 gs -1 .The light cylinder radius is defined as the radius 
where the corotational velocity is equal to the speed of light 

R\ c = —^ = 6.5 x 10 9 cm, ( 8 ) 

27r 

where in the last step the period of the ultra-luminous X-ray 
pulsar is substituted. When the Alfven radius is less than the 
light cylinder radius, the accretion flow may interact with 
the central neutron star. The rcorotation radius is defined as 
the radius where the corotational velocity equals the local 
Keplerian velocity 

/ GMP 2 \ 1/3 „ 1/3 

Rco = f 4n2 J = 1-8 x lO 8 !^ 73 cm, (9) 

the period of the ultra-luminous X-ray pulsar is substituted. 
When the Alfven radius is less than the corotation radius, 
large scale accretion is possible. Otherwise, the accretion 
flow is stopped by the centrifugal barrier (i.e. propeller, 
Illarionov & Sunyaev 1975). Accretion equilibrium corre¬ 
sponds to when the corotational radius is equal to the Alfven 
radius: R co = Ra (Lai 2014). The dipole magnetic field can 
be inferred from the spin-up trend of the ultra-luminous X- 
ray pulsar (Tong 2015c). The angular momentum carried by 
the accreted matter onto the neutron star is M acc y/GMR.A■ 
The spin-up rate of the central neutron star is determined by 

- 2 ttI-^ = MsccVGMRa, ( 10 ) 

where I is the star’s moment of inertia. The mass accretion 
rate can be inferred from the X-ray luminosity considering 
possible beaming effect. For a two solar mass neutron star, 
with a beaming factor of 0 . 2 , the dipole magnetic field is 
about (Tong 2015c) 

B P = L 6 x 10 12 i? 6 Lr o 3 40 G, (11) 

where R§ is the neutron star radius in units of 10 6 cm, 
-Liso, 4 o is the isotropic luminosity in units of 10 40 ergs -1 . 

Several other interesting aspects can also be obtained 
(Tong 2015c): (1) M82 X-2 is a transient because it may 
switch between accretion phase and propeller phase. ( 2 ) 
The soft X-ray excess in M82 X-2 may due to emission 
from the accretion disk. Future more accurate determina¬ 
tion of the inner disk radius will constrain different models. 
(3) The theoretical pulsation period of accreting magnetars 
may be very wide which ranges from 0.1 s to 10 2 3 * s. (4) Up 
to now, three signatures of accreting magnetar are available: 
magnetar-like bursts, a hard X-ray tail, and an ultraluminous 
X-ray pulsar. 


5 Discussion 

Both the radiation and rotation behaviors of NuSTAR 
J095551+6940.8 may due to effects of magnetic field. Ac¬ 
cording to the current models, there are at least three ways 
to constrain the magnetic field strengt!@: 

1. Maximum luminosity in the presence of a strong mag¬ 
netic field. The observed isotropic X-ray luminosity 
L x ft: 10 4 o ergs -1 should be smaller than the maxi¬ 
mum accretion luminosity onto magnetic neutron stars. 
The maximum luminosity depends on the magnetic 
field strength and has been modeled by several authors 
(Paczynski 1992; Mushtukov et al. 2015). 

2. Equilibrium state. The X-ray luminosity may switch 
from accretion phase to propeller phase (Tong 2015c; 
Dall’Osso et al. 2015; Tsygankov et al. 2015). The pul¬ 
sar may be in spin equilibrium (Eksi et al. 2015). Both 
of these two cases correspond to accretion equilibrium 
state which is characteristized by equaling the corota¬ 
tion radius and Alfven radius. 

3. Modeling the rotational evolution. The spin-up trend 
of NuSTAR J095551+6940.8 depends on the accretion 
torque. By assuming some kind accretion torque, the 
magnetic field strength can be obtained (Bachetti et al. 
2014; Eksi et al. 2015; Tong 2015c; Dall’Osso et al. 
2015). 

During the calculations, there are various uncertainties: the 
exact value of the X-ray luminosity (it changes with time); 
the beaming factor; the form of accretion torque; and vari¬ 
ous factors of 0.5 etc. Therefore, different authors may give 
different results of magnetic field strength. The exact value 
of magnetic field strength may range from 10 12 Gtol 0 14 G. 
It seems that the behavior of NuSTAR J095551+6940.8 can 
be understood in a very wide parameter space at present. 
Future more observations of more sources will make clear 
the theoretical uncertainties. 

The consequences of an ultra-luminous X-ray pulsar 
phase will be significant, especially for the star mass, mag¬ 
netic field, and rotation period. For an ultra-luminous X-ray 
pulsar phase with luminosity L x ~ 10 4 O ergs -1 , the cor¬ 
responding mass accretion rate will be M ~ 10 2 O gs -1 . 
For a neutron star with initial mass of 1.4 Mg, it will grow 
to 2 Mq in about 4 x 10 5 yr. Its mass will be 2.5 M 0 in 
about 7 x 10 5 yr. Therefore, the ultra-luminous X-ray pulsar 
phase is very efficient in producing massive neutron stars if 
this phase can lasts about 10 6 yr (Fragos et al. 2015; Wik- 
torowicz et al. 2015; Guo et al. 2015bIn analogy with 
classical accreting neutron stars, the star’s magnetic field 
may be reduced by the accretion matter (equation ( 20 ) in 
Zhang & Kojima 2006 and references therein). For an ini¬ 
tial magnetar strength dipole field of 10 14 G, after accreting 
about 0.6 Mq, the final magnetic field will be about 10 9 G- 
10 11 G. Assuming spin equilibrium, the star’s final rotation 

2 Measurement of cyclotron lines in the future will be very useful. 

3 The massive pulsar itself may account of the super-Eddington lumi¬ 

nosity of the ultra-luminous X-ray pulsar (Guo et al. 2015b) 
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period will be about 0.3 ms-14 ms (not considering vari¬ 
ous instabilities which may limit the minimum spin period). 
Therefore, the ultra-luminous X-ray pulsar phase may be 
another channel to form millisecond pulsarO (Pan et al. in 
preparation). If sub-millisecond pulsars exist, they can be 
formed via the ultra-luminous X-ray pulsar phase. 

6 Conclusion 

From the studies of isolated magnetars, it is known that the 
difference between normal pulsars and magnetars is their 
multipole field. A neutron star with a strong dipole field 
only is not a magnetar. Super-slow X-ray pulsars are prob¬ 
ably not accreting magnetars. They may be accreting high 
magnetic field neutron stars. The ultra-luminous X-ray pul¬ 
sar NuSTAR J095551+6940.8 may be an accreting magne¬ 
tar. Its radiation and rotation behaviors can be understood 
in a wide parameter space. However, the exact value of its 
magnetic field strength is still uncertain. It may be an ac¬ 
creting low magnetic field magnetar with multipole field 
of 10 14 G and dipole field of 10 12 G. This point of view 
is consistent with the study of isolated magnetars which 
tells that an aged magnetar is more likely to be a low mag¬ 
netic field magnetar. An ultra-luminous X-ray pulsar phase 
in the binary evolution may result in massive millisecond 
pulsars. Three signature of accreting magnetars are known 
at present: (1) magnetar-like bursts; (2) a hard X-ray tail; 
and (3) an ultra-luminous X-ray pulsar. Future hard X-ray 
telescopes (e.g. the HXMT telescope in China) may con¬ 
tribute to these aspects. 
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